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We present a family of exactly solvable spin- A quantum hamiltonians on a 3D lattice. The degenerate ground 
state of the system is characterized by a quantum error correcting code whose number of encoded qubits are 
equal to the second Betti number of the manifold. These models 1) have solely local interactions, 2) admit a 
strong-weak duality relation with an Ising model on a dual lattice 3) have topological order in the ground state, 
some of which survive at finite temperature, 4) behave as classical memory at finite temperature. The associated 
quantum error correcting codes are all non-CSS stabilizer codes. 
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I. INTRODUCTION 

One of the motivations for studying quantum error correct- 
ing code on lattice is to protect quantum information with- 
out active correction. Many models on 2D lattices have been 
proposed and analyzed — — but no-go theorem rules out all 
finite-range finite-strength hamiltonian system in 2D as a self- 
correcting quantum memory. 8 - 9 This does not apply to higher 
dimensions. For instance, it was shown that 4D toric code is a 
self-correcting quantum memory ! 10,11 Bombin et al. showed 
that there is also a 6D model that exhibits similar behavior. 12 
Whether such thermally protected model exists in 3D remains 
as an open problem. 3D toric code can store classical infor- 
mation at finite temperature but it fails to do so for quan- 
tum information. 13 Toplogical color code in 3D, albeit lack- 
ing a rigorous proof, is believed to show a similar behavior: 
there exists a string-like logical operator which is thermally 
unstable. 14 3D model proposed by Nussinov and Ortiz shows 
similar behavior. Another model was proposed by Cha- 
mon and analyzed recently by Bravyi et al. This model may 
be able to protect quantum information, but not in a thermo- 
dynamic sense i 17 ! 18 

It is worth noting that all the listed 3D models except Cha- 
mon's model share a similar property: the quantum error cor- 
recting code defining the ground state of the system is a CSS 
code, meaning that it can be decomposed into two classical 
codes. When studying the stability of these models, one can 
show that one of the codes can protect classical information 
from thermal fluctuation while the other one cannot. This 
means that there is a manifest difference between how the 
models treat the bit flip error and the phase flip error. Cha- 
mon's model treats X, F, and Z error in an identical manner 
but it lacks stability in thermal sense. Since we expect a singu- 
lar behavior at the phase boundary between an 'ordered state' 
and 'disordered state' for thermally stable quantum memory, 
absence of finite-temperature phase transition seems trouble- 
some unless there is an argument that can evade this logic. 
Motivated by these ideas, we present a new spin-^ model 
with finite temperature phase transition whose ground state 
is a non-CSS quantum error correcting code. Our model ex- 
hibits a topological order, but only the classical part survives 
in finite temperature. 

The outline of the paper is as follows. We set the stage by 
introducing the hamiltonian in Section HI] In Section [HI] we 
study the quantum code that defines the ground state of the 



hamiltonian. We calculate the number of qubits and find the 
logical operators. In Section lTVl we study the low-energy ex- 
citation of the hamiltonian that consists of particles and closed 
strings. We construct a duality relation with classical Ising 
model in Section[V]to show the finite temperature phase tran- 
sition. 



n. MODEL 

We place qubits on a vertices of a 4-valent 3D lattice. Using 
the notation Xi = a*, F, = crj, Zj = o; stabilizer generators are 
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where p is the plaquette and {i E p} denotes a set of ver- 
tices on plaquette p. We shall partition a set of plaquettes 
into P x ,Py,P z , which corresponds to a set of nontrivial sup- 
ports for B x p ,B y Pl B z p . We shall call elements of these sets as 
X— , Y—,Z— plaquettes. 

Our model is inspired by the construction of topological 
color code in 3D,™ For this quantum code, qubits reside on 
the vertices of the lattice, and the lattice is locally 4-valent. 
The stabilizer generators are either a product of Xs or product 
of Zs, and they correspond to the unit cells of different dimen- 
sions; in one example, generators are either in cubic form or 
plaquette form. Our approach differs in a sense that we only 
allow plaquette operators as stabilizer generators. 




(a) Vertex Figure 



(b) Unit Cell 



FIG. 1: Vertex figure and unit cell of our model. Qubits reside on the 
vertices. One can see that meets with another B x p ^ at one vertex 
whereas it meets with B y Py and B z p at two vertices. 
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Local description of our model can be seen in FIG |l(a)| At 
each vertex, there are 6 plaquette operators that have nontriv- 
ial support on it. Each plaquette operators meet with a same 
kind of plaquette operator on each vertices and meet with 4 
other plaquette operators on 2 vertices. Thus the assignment 
in FIG | 1 (a)| guarantees commutativity between the stabilizer 
operators. We must point out that not every lattice structure al- 
lows vertex figure like FIG |l(a)| There are only 4 translation- 
ally invariant convex tessellations that have tetrahedral ver- 
tex figure: bitruncated qubic honeycomb, cantitruncated cubic 
honeycomb, omnitruncated cubic honeycomb, and cantitrun- 
cated alternated cubic honeycomb. 19 Only the first three ad- 
mits an arrangement of plaquette operators similar to FIG |l(a)| 
at every vertex. In this paper, we mainly study the bitrun- 
cated qubic honeycomb model for its simplicity but analo- 
gous results shall be discussed in full generality if possible. 
Unit cell is shown in FIG |l(b)| and tessellation is shown in 
FIG|2] Bitruncated qubic honeycomb is a space-filling tes- 
sellation made up of truncated octahedra. It has 14 faces, 36 
edges, and 24 vertices. There are 6 square faces and 8 hexag- 
onal faces. Without loss of generality, one can set the 6 square 
faces to be Y plaquette operator, 4 of the hexagonal faces to 
be X plaquette operator and 4 remaining hexagonal faces to be 
Z plaquette operators. Hamiltonian is a sum over the plaquette 
operators. 

H = -J( £ B; x + £ B> py + £ B z pz ), (2.4) 

p x €P x p y eP y p z eP z 
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FIG. 2: Arrangement of stabilizer generators. Translation of unit 
cells form a tessellation. 



III. QUANTUM CODE 

Purpose of this section is to study the quantum code gener- 
ated by a set of group generators {B x p ^,B y Py ,B~ p _}. The section 
is mainly divided into two parts. In Section Itl! Al we count 
the number of encoded qubits. In Section lTHBI we completely 
specify a set of logical operators for each qubits. 



A. Number of Encoded Qubits 

Number of encoded qubits can be computed from the size 
of the stabilizer group and the number of physical qubits. 
Since the plaquette operators are not independent to each 
other, we must count the number of independent relations. 
In such pursuit, geoemetrical interpretation of our model be- 
comes useful. We would first like to point out that multiplying 
all the plaquette operators on a unit cell reduces to identity. 
One can see this from FIG |l(b)| Since any contractible closed 
surface on the lattice can be represented as a union of unit 
cells, one can see that multiplication of plaquette operators on 
any contractible closed surface reduces to identity. Therefore 
we have C — 1 independent relations which generate smooth 
deformation, where C is the number of 3-cells. We must sub- 
tract 1 becuase multiplying all but one cell results in a relation 
for that very cell. 

Let us consider a periodic boundary condition on all 3 di- 
rections. There exists noncontractible surface that reduces to 
identity as one can see in FIG |3(a)| FIG |3(b)| Since there are 3 
topologically distinct noncontractible surfaces, we have 3 in- 
dependent relations, resulting in C + 2 independent relations. 
Finally, multiplying all X-like operators adds one independent 
relation. One can check that multiplication of Ys and multipli- 
cation of Zs are implied by the previously mentioned relations. 




(a) Top View (b) Side View 



FIG. 3: Representation of nontrivial constraints between the stabi- 
lizer operators. One can see that multiplication of all the plaquette 
operators on a noncontractible closed surface reduces to identity. At 
each vertex, there are either 1) exactly one X, one Y, and one Z or 2) 
two Xs and two Zs. 

Accounting for these relations, number of encoded qubits 
is V — F + C + 3 = 3. This reasoning can be generalized to 
any orientable 3-manifold. 

Lemma 1 For stabilizer group {B x p ^,B } Py ,B p }, k = £>2- 

Proof : We use the definition of Euler Characteristic. 

% = V -E + F-C = (3.1) 

% is trivially due to Poincare Duality. In the dual lattice, V 
is the number of tetrahedral cells. E is number of faces, and 
hence E = 2V . Therefore we have 

V-(F-C)=0. (3.2) 
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Hence 

k = V-(F-(C-l + l+b 2 )) (3.3) 
= b 2 , (3.4) 

where Z?2 is the second Betti number of the manifold. One can 
also use this intuition to prove that the group generated by the 
plaquette operators does not contain —I. 

Lemma 2 (B~ Px ,B y Pyl B z p _) does not contain —I. 

Proof: Any constraint between the plaquette operators can be 
represented as a product of closed 3-manifold. For each unit 
cell, we have 24 vertices at which X, Y, and Z meets. Since all 
the generators commute with each other, we can arrange the 
product to be the following canonical form. 

II„ «; ; . (3.5) 

Since XYZ = i, the product of plaquette operators on a unit 
cell is 1 . Similarly, product of plaquette operators on a non- 
contractible surface described in FIG |3(a)| FIG |3(b)( we have 
An vertices where X,Y, and Z meets. Hence we arrive at the 
same conclusion. Since any product of plaquette operators 
that results in a trivial operator can be constructed by these 
constraints, the group does not contain —I. 



B. Logical Operators 

There are two logical operators that are reminiscent to the 
surface and string operator of 3D toric code. These are drawn 
in FIG|U One can see the surface operator on the top of the 
lattice system which is a product of s on one layer of Y- 
plaquettes. The complementary logical operator to this is the 
string operator that has a sequence of YZYXYZYXYZYX ■ ■ ■ 
along the line perpendicular to the surface operator. This 
string winds around the torus and completes a noncontractible 
loop. These two operators anticommute with each other and 
both of them commute with the stabilizer generators. We can 
similarly define two sets of complementary operators in other 
directions. One can easily check the expected commutation 
and anticommutation relations. 



IV. LOW ENERGY EXCITATION 

Quasiparticles excitations in 2D typically arise as anyons. 
For instance, in Kitaev's toric code, two quasiparticles are cre- 
ated in pair, and when fused together, they vanish^- There 
are two kind of particles analogous to electric and magnetic 
charge, and when one particle winds around another one, the 
system attains a nontrivial global phase. In 3D, trajectory of 
winding around another particle can be deformed into a trivial 
contour. Hence one needs higher dimensional object to attain 
a similar topological action. In 3D there are closed string-like 
excitations and particle-like excitations. 2 ' 13 When the particle 
winds around the string so that the trajectory and the string 




FIG. 4: There is one surface operator and one string operator for each 
qubits. Surface operator corresponds to the product of ZZZZ on Y- 
plaquettes. String operator is the line perpendicular to this surface, 
showing a sequence YZYXYZYX 



together forms a knot, the system attains a nontrivial global 
phase. 

Our model presents a similar picture. Particle-like excita- 
tions are created in pair. If we truncate a string-like logical op- 
erator, excitations form at the end points. When the particle- 
antiparticle pair is created, they can diffuse without any extra 
energy cost. Closed string-like excitations can be similarly 
thought as a truncated surface-like logical operator. Near the 
boundary of the surface, there are excitations and hence the 
energy cost grows linearly with the size of the surface. When 
a particle penetrates the closed string, we find that 




FIG. 5: Representation of particle penetrating through a string-like 
excitation. Truncated surface operator is a product of Z-plaquettes in 
white. Trajectory of the particle is a nontrivial support of the colored 
plaquette operators, which coincides with the Z-surface. 
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tyb*ua)=SP\*) (4.1) 
WFinal) =USP\G>) = - \\fi nma l ) , (4.2) 

where S is a closed-string excitation, Pisa particle excitation, 
and U is a trajectory of the particle. Thus system gains e ,% 
phase factor. This is illustrated in FIG[5] One can see that as 
a particle penetrates through the surface operator and returns 
to the original position, it coincides with the surface operator 
at one vertex, thus giving the anticommutation relation. 

Low energy excitation in terms of elementary objects pro- 
vides us an intuitive picture for the thermal stability. Particles 
can be created out of vacuum in pair and propagate freely. 
They can diffuse and wind around the torus to induce logical 
error. Closed strings, on the other hand, need energy that is 
proportional to its perimeter. Given a closed string-like ex- 
citation as in FIGJU the stabilizer generators anticommuting 
with the surface operator only reside near the boundary of the 
surface. Z-plaquettes trivially commute with the surface op- 
erator. X -plaquettes commute with the surface operator since 
they meet at two vertices. However, there are T -plaquettes 
meeting at exactly one vertex at the boundary. Hence we ex- 
pect our system to be a stable classical memory. 

V. DUALITY 

Typical strong-weak duality relation relates a strong cou- 
pling limit of one model to a weak coupling limit of another 



£ c is a sum over a configuration of closed 2-manifolds. A c 
is the number of plaquettes for each configurations. C.T. cor- 
responds to the cross terms between closed 2-manifolds and 
space-filling product of Xs, Ys, or Zs. n Xi y >z corresponds to 
the number of X,Y,Z— plaquette operators. The main idea is 
that the partition function is dominated by the first term in the 
thermodynamic limit. We show this in AppendixlAl 

Lemma3 Z-C.T. - (oc % +OC n ? +0C n *) =Z IG $J), whereZ 1G 
is a partition function oflsing gauge theory on the same lattice 
with temperature p and coupling constant J. 

Proof : Consider a mapping B* — > ZZZZZZ, B y Py — s- ZZZZ, 
B z p _ — > ZZZZZZ, where Z ■ • • Z are products of Z on the edges 
of each plaquettes. The resulting model is an Ising gauge the- 
ory on a bitrucated cubic honeycomb. Partition function is 

Z /G =fr(exp(-p//)) (5.6) 
= (coshp7)"rr(l +tanhp/5,-), (5.7) 



model: we use a slightly different strategy here. We first show 
that our model can be mapped into an Ising gauge theory, from 
which we can use the Wegner-type duality relation with Ising 
model. Mapping from our model to Ising gauge theory is not 
exact for finite sized lattice, but this difference vanishes in the 
thermodynamic limit. Starting from the partition function of 
our model, 

Z = fr(«p(-Pff)) (5.1) 
= t r(n Ste s (cosh p7 + 5,- sinh p/)) , (5.2) 

where5,G{B; t ,4 v ,^}, 



Z=(coshp/)V(IIi(l+aSj)) (5.3) 
l 

= (coshpj)V( £ iW'Sf)- (5-4) 

M=o 

Since the Pauli operators are traceless, the nonvanish- 
ing terms correspond to the nontrivial constraints presented 
in Section IIII Al Note that there were two kind of con- 
straints: constraints coming from the closed 2-manifold and 
constraints coming from space-filling products of X, Ts, or 
Zs. Using this, we can write down the partition function in the 
following form. 



(5.5) 



where 5,-s are either ZZZZZZ or ZZZZ depending on the pla- 
quette. Since Pauli operators are traceless, only a product of 
plaquette operators that are union of closed surface survives. 
Therefore, we conclude 

Z /G (P7) = Z C.T. - (a"> + a"y + a"=). (5.8) 

Using the duality relation between Ising gauge theory and 
Ising model, we can map our model into an Ising model. We 
show the duality relation in AppendixlBl 

Theorem 1 Our model with coupling constant fij is dual to 
the classical Ising model on a dual lattice with a dual coupling 
constant $J = — j lntanhp/. 

Since the Ising model undergoes a finite temperature phase 
transition, so does our model. This is analogous to the be- 
havior of 3D toric code under temperature change. As in our 
model, one can show that 3D toric code has critical temper- 
ature by using the duality relation with Ising model. Below 
the critical temperature, there is a symmetry breaking with re- 
spect to a surface-like logical operator. Symmetry associated 



J 



Z = (2coshpj)"(£a Ae + (1 +a"- r )(l + a">')(l + a" : ) - 1 +C.T.) 

C 

I 
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to the string-like logical operator is broken only at the ground 
state. 

One glaring difference though, is that 3D toric code can be 
decomposed into two classical hamiltonians without spoiling 
the phase transition: the hamiltonian responsible for correct- 
ing the bit flip error is identical to Ising gauge theory, which 
has finite temperature phase transition. On the other hand, 
the hamiltonian responsible for correcting the phase flip error 
does not have a phase transition. Hence one can intuitively un- 
derstand that 3D toric code can only correct bit flip errors but 
not phase flip errors under thermal equilibrium. Our model 
does not allow such decomposition. Once we get rid of any of 
Bp , B Py , or B z p , the partition function does not exhibit a phase 
transition any more. This shows that non-CSS code with finite 
temperature phase transition in 3D does not necessarily pro- 
vide a self-correcting quantum memory. 



VI. CONCLUSION 

In this paper, we studied an exactly solvable 3D spin model 
and studied its topological order. The ground state of the sys- 
tem defines a non-CSS quantum error correcting code. At fi- 
nite temperature, this system is expected to behave as a stable 
classical memory, but not as a stable quantum memory. This 
is mainly due to the fact that there exists a string-like logical 
operator. In light of studying the possibility of self-correcting 
quantum memory, this reconfirms the general properties that 
have been found in 3D stabilizer codes so far: for each en- 
coded qubit, there exists one surface-like logical operator and 
one string-like logical operator. It seems that we cannot avoid 
such outcome unless the shape of the logical operator changes 
as the system size changes, as in Chamon's mode L 17 ' 18 This 
in fact was recently argued to be the general feature of stabi- 
lizer codes whose number of encoded qubits remain invariant 
under system size changed 

It is worth noting that the thermal stability analysis of our 
model is not rigorous at this stage, even though the energy bar- 
rier increasing as the perimeter of the surface is a compelling 
evidence that this must be true. It would be desirable to make a 
rigorous estimate of thermal relaxation rate using the method 
introduced by Chesi et al. 21 We expect the string-like logical 
operator to be thermally fragile and the surface-like logical 
operator to be stable. As in 3D toric coder^ we also expect 
the topological entropy of our model to show a singular behav- 
ior near the critical point. These singular behavior arise due 
to the existence of finite temperature phase transition, which 
we can show rigorously by the strong-weak duality relation 
between our quantum model to a classical Ising model on the 
dual lattice. 
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Appendix A: Bound for the cross terms. 

Corss term can be written as 



a' 



fii—lnS 



(Al) 



i£{x,y,z} 



where n x ,n y ,n z are total number of X,Y,Z— plaquettes and 
n x ,n y) n z are number of X,Y,Z— plaquettes for configuration 
c. 

Lemma 4 There exists < £1,2 < 1 such that 

A c + rij - 2n c ( > Z\A C + £ 2 rii (A2) 

for Vc, i. 

Proof : Consider i = x. Left hand side of the inequality is 

n° y + n\ -n L x + n x > n° y + n\ - (1 - e)n c x + (1 - e)n x (A3) 



>{^)A c + (l-e)n x 



(A4) 



On the second line, we used the fact that the minimum is 
achieved in the case where n y = 0, implying n\' = n x = \A C . 
Same logic can be applied to i = z. For i — y, 

n c x + n\ - n c x + n y >n x +n c z -(l- e)n c y + (1 - e)n y (A5) 

>(~-~(l-e))A c + (l-e) ra> . (A6) 

Similarly, here we used the fact that the minimum is achieved 
in the case where one of n c x or n c z is 0. Then we have a 2 : 3 ra- 
tio between the X — (Z— )plaquettes and Y — plaquettes. There- 
fore, for e > A, we have such (£1,62). 



Lemma 5 



hm — -tk^t -> 1- 
m l^~Z IG {$J) 



(Al) 



, where Z/c(p/) is a partition function for Ising gauge theory 
with temperature (3 and coupling constant J. vol is the volume 
of the lattice. 



Proof : 

We use 



L f , (2 cosh (3/')" ^ 1 a /Aos 
a ZlAc = ) ) a' Ac (A8) 

(2coshp7')"'T 



1 



2 cosh (3/ 



7 )"Z/g(PA 



where 



tanhp7' = (tanhp7) £l . 



(A9) 



(A10) 
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Thus the cross terms can be bound by 



Proof: 



(All) 



where 8, = ^,where n is the total number of plaquettes. This 
becomes 



z /G (p/)((^4-)^r, 



i -r E i 



(A12) 



_ 2 1 £ 2 

where f = tanhp/'. One can show that ( t 1 )?t I i B < 1 for 

i-t*i 

p7 > 0. Since the renormalized coupling constant /' is larger 
than J, we can see that these correction terms become negli- 
gible in thermodyamic limit. Therefore, 

where J' > J and < A. < 1. In n — > °° limit, we get the desired 
result. 



Appendix B: Duality between Ising gauge theory and Ising 
model 



Lemma 6 Ising gauge theory on bitruncated cubic honey- 
comb is dual to Ising model on the dual lattice. 



Z = (coshpy)"fr(n,(l +tanhp/S 1 -)) (Bl) 
l 

= (coshp/)"fr( £ iW'Sf') (B2) 
l 

= (2coshp/)" £ rLoc*<n e 8 2 (£> ;> ), (B3) 

ft}=o i 

where Yl e is a product over all the edges and Y*jky,e is a 
sum over kjs that have nontrivial support on edge e. There 
are three such kjS. One can use kj ;e = 3(1 —ZZ), where 
ZZ is a product of Zs on qubits that reside on the vertices 
of the dual lattice. For 8 spin configurations (Zi.Z^.Z^) = 
(-1,-1,-1), (1,1,1), (1,-1,-1), (-1,1,-1), (-1,-1,1), 
(1,1,-1), (—1,1,1), (1,-1, 1), one can see that all of these 
configurations satisfy the delta function. Furthermore, we 
have 2 combinations for (k\,k2,kj,) = (0,0,0), 2 combinations 
for (0,1,1), (1,0,1), and (1,1,0). Plugging this in, we get 



1 , 

z=(coshpy)" £ nia}~^ Zi +"i Zi -"i, 



(B4) 



where Z;±/j ; is the Z operator on the dual sites of plaquette i. 
hi is the unit normal vector to the plaquette. Therefore, up to a 
constant, partition function is identical to the partition of Ising 
model with p/ = — \ lntanhp/. 
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